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: Nicotine addiction is associated with risky behaviors and abnormalities in local brain areas related 

: to risky decision-making such as the dorsal anterior cingulate cortex (dACC), anterior insula (Al), and 

: thalamus. Although these brain abnormalities are anatomically separated, they may in fact belong to 

: one neural network. However, it is unclear whether circuit-level abnormalities lead to risky decision- 
making in smokers. In the current study, we used task-based functional magnetic resonance imaging 

: (fMRI) and examined resting-state functional connectivity (RSFC) to study how connectivity between 

: the dACC, insula, and thalamus influence risky decision-making in nicotine addicts. We found that 

: an increase in risky decision-making was associated with stronger nicotine dependence and stronger 
RSFC of the dACC-rAl (right Al), the dACC-thalamus, the dACC-IAI (left Al), and the rAI-IAI, but that 
risky decision-making was not associated with risk level-related activation. Furthermore, the severity 

: of nicotine dependence positively correlated with RSFC of the dACC-thalamus but was not associated 

: with risk level-related activation. Importantly, the dACC-thalamus coupling fully mediated the effect 
of nicotine-dependent severity on risky decision-making. These results suggest that circuit-level 
connectivity may be a critical neural link between risky decision-making and severity of nicotine 
dependence in smokers. 


: Nicotine addiction is associated with risky behaviors!” and abnormalities in brain areas related to decision mak- 
: ing* such as the dorsal anterior cingulate cortex (dACC)**, insula”*, thalamus”? and striatum!)””, 

The dACC is primarily responsible for evaluation in decision making”. Abnormal activity in the dACC 
: has been linked to cravings for cigarettes’. The anterior insula (AI) plays a key role in mediating the relation- 
: ship between nicotine dependence and risky decision making [e.g., Risky-Gains task'*]. Abnormal activity in 
: the thalamus has also been associated with cigarette cravings’’. Although brain abnormalities associated with 
: risky decision making and nicotine addiction are found in anatomically separated regions (e.g., the dACC, AI 
: and thalamus), these regions may in fact be parts of one neural network. Changes in this neural network may 
: reveal important information about global deficits in the organization of the brain'®. So far, it is unclear whether 
: circuit-level abnormalities lead to behavioral deficit in risky decision making in population of nicotine addicts. 
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Figure 2. The mean pumps positively correlated with nicotine dependence severity. 


The dACC, AI and thalamus are critical nodes of the salience network, a well-defined functional brain net- 
work” that mainly serves to identify relevant internal and external stimuli to guide behavior". Risky stimuli (e.g., 
a risky choice), which are associated with strong subjective excitement and strong interoceptive change”, may be 
highly salient. Previous studies found stronger connectivity within the salience network in nicotine addiction”). 
In the present study, we hypothesized that in population of nicotine addicts, greater connectivity between the 
dACC, AI, and thalamus may enhance the salience processing of risky stimuli, resulting in an increase in risky 
decision making. 

To study how functional connectivity influence risky decision making in population of nicotine addicts, we 
used task-based functional magnetic resonance imaging (fMRI) in conjunction with the Balloon Analogue Risk 
Task (BART Fig. 1) and resting-state functional connectivity (RSFC), which can acquire temporal correlations 
between spontaneous regional activities in the absence of externally controlled tasks or stimuli”. We first tested 
the relationship between the severity of nicotine dependence and risky decision making using BART. Then, we 
used correlation and causal mediation analyses, statistical methods that can inform understanding of how brain 
regions interact to result in behavior, thus, whether the RSFC was the link between risky decision making and the 
severity of nicotine dependence. 


Results 

Task performance. The BART used in this study is a sequential decision-making task, widely used as a tool 
for assessing risky behavior”. It provides participants a chance to pump virtual balloons to win potential mone- 
tary reward. However, each pump simultaneously increases the risk to lose accumulated reward for that balloon”. 
The mean pumps (the mean number of pumps across trials for balloons that did not explode, which is a typical 
behavior index of risky decision making in the BART) were 5.7 + 1.4 and the FIND (Fagerstrom Test of Nicotine 
Dependence, which is used to measure the severity of nicotine dependence) score was 4.2 + 1.9. The mean pumps 
were positively correlated with the FIND scores across participants (r = 0.442, p = 0.045; Fig. 2). 


Relationship of task-based activation to task performance and severity of nicotine depend- 
ence. During pumping, significant risk level-related activations (p < 0.005, family-wise error corrected) were 
found in the bilateral dorsal anterior cingulate cortex (dACC), right anterior insula (rAI), left anterior insula 
(AT), bilateral thalamus, right precentral gyrus, right inferior parietal lobe (rIPL) and left inferior parietal lobe 
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| Cluster | Meanpumps | FTND 
Regions x y | z | Maxz| C e | p |r Pp 
R precentral gyrus —34 14 52 5.35 249 0.288 | 0.205 | —0.008 | 0.972 
R anterior insula —32 | —23 =l 6.65 146 0.399 | 0.074 0.192 0.405 
Bilateral dACC —8 —2 47 7.20 132 0.368 | 0.101 0.180 0.435 
RIPL —56 25 32 4.59 82 0.244 | 0.286 0.034 0.883 
Bilateral thalamus —8 22 —4 | 3.60 63 0.305 | 0.179 0.152 0.511 
Lanterior insula 31 —17 11 6.07 42 0.335 | 0.138 0.260 0.255 
LIPL 55 28 35 4.47 40 0.257 | 0.262 0.011 0.961 


Table 1. The risk level-related activation derived from a whole-brain analysis and its relationship to the 
mean pumps and severity of nicotine dependence. 


(IIPL) (Table 1). None of the activations of these regions significantly correlated with the mean pumps or FIND 
scores (Table 1). 


RSFC and relationship to task performance and severity of nicotine dependence. The RSFCs 
between the bilateral dACC, r-AI, 1-AI and bilateral thalamus (Fig. 3) were calculated and the relationship to 
the mean pumps and FTND score was tested. Mean pumps positively correlated with the RSFC of dACC-rAI 
(r = 0.525, p= 0.015), dACC-thalamus (r = 0.576, p = 0.006), dACC-IAI (r = 0.564, p = 0.008), and rAI-lAI 
(r = 0.465, p = 0.034), but did not correlate with the other RSFC. The FIND score positively correlated with 
the RSFC of the dACC-rAI (r = 0.473, p = 0.031) and dACC-thalamus (r = 0.647, p = 0.002), but did not corre- 
late with the other RSFC. We performed a FDR-correction for multiple comparisons (FDR < 0.05). Correlation 
between the mean pumps and the RSFC of dACC-thalamus (r = 0.576, p = 0.036, corrected), dACC-1AI 
(r=0.564, p = 0.024, corrected), dACC-rAI (r = 0.525, p = 0.030, corrected), and rAI-lAI (r = 0.465, p= 0.051, 
corrected), and correlation between the FTND score and the RSFC of dACC-thalamus (r = 0.647, p = 0.012, 
corrected) survived the correction. However, correlation between the FTND score and the RSFC of dACC-rAI 
(r= 0.473, p = 0.093, corrected) did not survive. Only the RSFC of the dACC-thalamus correlated with both the 
mean pumps and the FIND score (Fig. 4). 

Other RSFC between regions related to risk level except couplings between bilateral dACC, bilateral AI and 
bilateral thalamus were calculated and then relationship to mean pumps and FIND score was tested (see details 
in supplementary materials). We did not find any coupling that significantly correlated with both mean pumps 
and FIND score. 


RSFC mediation effect. Mediation analysis is a statistical method widely used in psychology research 
to examine correlated pathways among relational variables and tests whether association between an inde- 
pendent and dependent variable could be explained by a third variable’. To examine whether the RSFC of the 
dACC-thalamus pathway might mediate the effect of the FIND score on mean pumps, we performed a medi- 
ation analysis”*?> (see details in methods section). The mediation effect was considered to occur if the relation- 
ship between the FIND score (independent variable) and mean pump (dependent variable) could be partially 
or totally accounted for by the RSFC of the dACC-thalamus (hypothesized mediator). In the first step of the 
mediation analysis, the hypothesized mediator of the RSFC of dACC-thalamus was regressed on the FTND 
scores. Higher FIND score predicted a stronger RSFC of dACC-thalamus [8 = 0.647, SE = 0.015, t(19) = 3.694, 
p=0.001]. Secondly, the mean pumps were regressed on the FIND score, and the higher FIND score signifi- 
cantly predicted more mean pumps [8 = 0.442, SE = 0.275, t(19) = 2.145, p = 0.023]. Thirdly, the mean pumps 
were regressed on the FIND score together with the hypothesized mediator. Stronger RSFC of dACC-thalamus 
significantly predicted more mean pumps [8 = 0.499, SE = 2.151, t(18) = 1.990, p= 0.031]. The relationship 
between the FTND score and mean pumps decreased in strength with the presence of the mediator (RSFC of 
dACC-thalamus) and was not significant [8 = 0.119, SE = 0.188, t(18) = 0.473, p = 0.321]. The mediation analysis 
indicated that the RSFC of dACC-thalamus fully mediated the effect of the severity of nicotine dependence on 
the mean pumps (Fig. 5). A partial correlation analysis was employed to support the mediation effect. The FIND 
score was not correlated with the mean pumps while controlling for the influence of RSFC in the dACC-thalamus 
pathway (r (18) = 0.111, p = 0.642) (see the detail in supplementary materials). 


Discussion 

In the present study, we found that the FIND was positively associated with mean pumps. More importantly, 
greater RSFC among the dACC, AI, and thalamus was associated with higher mean pump values and more 
severely nicotine dependence in smokers. 

Our results that showed association between greater RSFC of dACC-rAI, dACC-thalamus, dACC-IAI and 
rAI-IAI and increased risky decision making (i.e., the mean pumps) suggested that circuit-level engagement 
of the salience network may be required in risky decision making. One study has shown that the activity in a 
large network of brain regions tasked with processing risk during BART could predict choices in future risky 
decisions”. The dACC, AI, and thalamus have been suggested as important nodes of the salience network that 
respond to the degree of subjective salience, whether emotional, cognitive, or homeostatic”””*. Thus, it is possible 
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Figure 3. Significant risk level-related activation during risky decision making (regions of interest for 
RSFC). Panel (a). During risky decision making, risk level-related activation was found in the bilateral dorsal 
anterior cingulate cortex (dACC). Panel (b). Bilateral thalamus. Panel (c). Right anterior insula (rAI). Panel (d). 
Left anterior insula (1AT). 


that facilitated communication within the salience network is represented to enhance the processing of salient 
stimuli”. 

Risk level-related activations of brain areas were not associated with mean pumps in the present study, which 
was consistent with previous studies using BART?+°°-°. It is possible that because decision making employs a 
complex set of cognitive functions, circuit-level engagement possesses more explanatory power than the function 
of a single region. Furthermore, recent evidence suggests that decision making is mediated by a brain network 
centered on the dACC and that elements of the network are interconnected!*. Functional networks, rather than 
single brain nodes, have amassed the ability to support complex functions such as decision making**. 

Studies of resting-state? and task-based** functional connectivity in population of nicotine addicts have 
shown that connectivity with the dACC is modulated as a function of severity of nicotine dependence. However, 
whether the link between risky decision making and severity of nicotine dependence can be explained by changes 
in functional connectivity in population of nicotine addicts has not been directly tested. The present study found 
that the dACC-thalamus coupling fully mediated the relationship between severity of nicotine dependence and 
mean pumps in population of nicotine addicts. The partial correlation data showed consistent result (see the par- 
tial correlation analysis in supplemental materials). These results suggest that changes in the RSFC of the salience 
network (i.e., the dACC-thalamus) may be the potential underlying neural mechanism of risky decision making 
and nicotine dependence in population of nicotine addicts. The thalamus carries feedback from the basal ganglia 
to the dACC*”. The dACC, which is part of the salience network and motor planning circuits** and is associated 
with increased risk-seeking decisions”, initiates appropriate control signals to regulate behavior”. Combined 
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Figure 4. RSFC and its relationship to the FTND score and the mean pumps. The RSFC of the dACC- 
thalamus positively correlated with the FTND score and the mean pumps. 
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Figure 5. RSFC mediation effect. The RSFC of the dACC-thalamus fully mediated the effect of severity of 
nicotine dependence on the mean pumps. Values presented are the standardized regression coefficients. The 

value in parentheses represents the coefficient for the direct (i.e., unmediated) path. *p < 0.05, **p < 0.01, one- 
tailed. 


with our results, these findings suggest that severely dependent smokers might engage in enhanced approach 
motivation when seeking motivationally salient stimuli, leading to an increase in risky behavior. 

Although we also tested other RSFCs, we did not find any couplings that significantly correlated with both 
the mean pumps and the FIND score (see details in the supplemental materials). Two dissociable networks are 
involved in a large-scale brain network called the task-positive network“: one, termed the salience network, is 
anchored in the dACC and anterior insula, whereas the second is composed of the parietal and lateral prefron- 
tal regions and is termed the executive control network’. In this study, we found RSFCs that were significantly 
correlated with both the mean pumps and the FIND score only between key nodes of the salience network (i.e., 
RSFCs of dACC-thalamus). We did not observe these results within the executive control network, suggesting 
that a specific circuit-level change in the salience network might be the underlying neural mechanism of risky 
decision making and nicotine dependence in population of nicotine addicts. 

In addition, the RSFC of the dACC-thalamus positively correlated with severity of nicotine dependence. 
According to cognitive control deficit theory“’, poor self-control ability of drug addicts may result in drug use. 
It is therefore expected that the activation of regions of the frontal cortex such as the DLPFC*, which plays a key 
role in cognitive control”, should be inversely correlated with severity of nicotine dependence. According to 
the incentive sensitization theory” and the interoception theory“, severely dependent smokers show increased 
activation of the insula and dACC as compared to less-dependent smokers, which may represent increased inter- 
oceptive processing and motor preparation that lead to motivational drive to smoke**“. In the present study, we 
found positive correlation between the RSFC of the dACC-thalamus and severity of nicotine dependence, which 
may support the incentive sensitization theory and the interoception theory. However, we did not find significant 
task-based activation of the frontal cortex**”. But, the frontal cortex had a trend to activate (see the details in 
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supplementary materials), making it difficult to explain the cognitive control deficit theory, and more evidence is 
needed to support this theory. 

In previous studies, smokers took more pumps than non-smokers on the BART!™®. Thus, we inferred that 
smokers were more likely to take risks than non-smokers. In the population of nicotine addicts, we found a pos- 
itive correlation between the FTND and the mean pumps, suggesting that severity of nicotine dependence may 
influence risk taking. 

This study had several limitations. Firstly, the lack of a control group of non-smokers made it hard to 
directly assess whether smokers actually took more risks and had stronger connectivity of salience network 
than non-smokers. Our results showed that severity of nicotine dependence could influence risk taking behav- 
ior by influencing connectivity of salience network in population of nicotine addicts. Our results suggest that 
circuit-level changes in the salience network might be the underlying neural mechanism of nicotine addic- 
tion and abnormal risky decision making. However, more studies are required to conclusively ascertain this 
claim. Secondly, some decision-making phases, such as anticipation and selection, cannot be isolated by BART. 
However, the present study focused on the influence of risk level on risky decision making rather than on specific 
decision-making phases. Thirdly, only one female participated in the present study, so we could not analyze the 
role of gender. Though gender was not considered in this study, it might be an important factor and merits further 
research in the future. Fourthly, the difficulty of decision under inflation, which increased with increasing risk 
levels, could not be isolated. Difficulty is an important factor that should be considered before making a decision. 
In addition, anticipated rewards and losses are also considered in evaluation of risk-taking. While risk-taking is 
a result of evaluation of many variables, this specific variable (e.g., difficulty) was not considered in this study. 

In the present study, we found that circuit-level connectivity between key nodes of the salience network, which 
impacts risky decision making, is associated with severity of nicotine dependence. These results suggest that 
in population of nicotine addicts, circuit-level changes in the salience network might be the underlying neural 
mechanism of nicotine dependence and risky decision making. Thus, interventions, such as transcranial direct 
current stimulation and transcranial magnetic stimulation, to balance the functional connectivity of the salience 
network may enhance the treatments for risky decision making and nicotine addiction. 


Materials and Methods 

Participants. Internet advertisements were used to recruit healthy daily smokers (>10 cigarettes per day 
for at least 1.5 years). Twenty-one smokers (mean age: 26.38 + 4.31, range 18-40 years, 1 female) participated 
in the study. All participants had normal or corrected-to-normal vision and were right-handed. Potential par- 
ticipants were excluded if they reported a diagnosis of a psychiatric disorder, the use of addictive drugs (except 
nicotine), a prior head injury, or any contraindications for participating in an MRI study (e.g., non-remova- 
ble metallic implants). The severity of nicotine dependence was assessed using the Fagerstrom Test of Nicotine 
Dependence [FTND*’]. The study was approved by the Research Ethics Committee of the University of Science 
and Technology of China, and written informed consent in agreement with the Declaration of Helsinki was 
obtained from all participants. The methods and procedures were carried out in accordance with the approved 
guidelines. 


Task and procedure. The implementation of BART was modified from a previous study*. During this task, 
participants pressed buttons to sequentially pump a computer-simulated balloon image. Pressing the right button 
led to two possible outcomes that were immediately fed back to the participants: a no-explosion balloon pump- 
ing event in which a larger balloon with an increased wager was displayed, or a loss event in which a picture of 
an exploding balloon and the text “You Lose!” was displayed immediately. If the balloon exploded, participants 
would lose the wager and the virtual monetary amount lost was subtracted from the cumulative earnings as a 
penalty. If the participants chose to discontinue pumping (pressing the left button), a cash-out event in which the 
text “You Win!” was displayed and the amount of the reward was added to the cumulative earnings. 

There was a jittered 2-4s interval prior to the beginning of each balloon and a jittered 1.5-2.5 s interval 
between two continuous pumps. Twelve pump responses were possible for a given balloon, with a parametric 
increase in the probability of explosion over successive pump responses (i.e., take a risk with 0% probability of 
explosion for ¥0.0; 2.1% for ¥0.05; 4.2% for ¥0.15; 6.3% for ¥0.25; 14.6% for ¥0.55; 23.9% for ¥0.95; 31.3% for 
¥1.45; 43.8% for ¥2.05; 56.3% for ¥2.75; 68.8% for ¥3.45; 79.2% for ¥4.25; and 89.6% for ¥5.15). The maximum 
number of pumps and the exact probability of explosion associated with a given inflation were not provided to 
participants. The risk level in the BART was directly and ecologically defined as the probability of explosion for 
each balloon. 


Behavior analysis. Consistent with previous studies***!, risky decision making in this task is typically meas- 
ured as the mean number of pumps, or “mean pumps’, across trials for balloons that did not explode. To inves- 
tigate whether the nicotine addicts’ risky decision making was specifically affected by the severity of nicotine 
dependence, we used Pearson correlation analysis to compute the correlation between the FIND score and the 
mean pumps. 


MRI data acquisition. Gradient echo-planar magnetic resonance imaging data were obtained during the 
whole task procedure on a 3T Siemens Magnetom Trio scanner (Siemens Medical Solutions, Erlangen, Germany) 
at the Anhui Provincial Hospital. Task stimuli were projected onto a display screen at the back of the magnet’s 
bore, and participants viewed the stimuli through a mirror. Button-press responses on an {MRI-compatible 
response box were made using the left and the right thumbs. A circularly polarized head coil was used for radio 
frequency transmission and reception, with foam padding to restrict head motion. 
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Resting-state {MRI data (240 volumes, 8 minutes) were acquired by asking participants to keep their eyes 
closed and to remain at rest. This was followed by another 8 minute functional scan corresponding to BART. 
Functional images were acquired with a T2*-weighted echo-planar imaging sequence (TE= 30 ms, TR= 2000 ms, 
FOV = 240 mm, Matrix = 64 x 64, flip angle = 85°) with 33 axial slices (no gaps, voxel size: 3.75 x 3.75 x 3.7mm?) 
covering the whole brain. In addition, high-resolution T1-weighted spin-echo images were also acquired. 


Image data preprocessing. The imaging data were processed using AFNI [Analysis of Functional 
Neuroimages*’] and MATLAB (version 7.6.0.324; Mathworks). Each participant’s raw data were corrected for 
temporal shifts between slices, corrected for motion, spatially smoothed with a Gaussian kernel (full width at half 
maximum = 8 mm), and temporally normalized (for each voxel, the signal of each volume was divided by the 
temporally averaged signal). 


Analysis of functional imaging data. Three types of events were included in the general linear model 
of the BART task: pumps, cash outs and balloon explosions. Each type of event included two regressors in the 
general linear model. For the pumps, which represented risk-taking in this task, two regressors were defined: 
(1) Pumps, arametrico the only regressor we were interested in, which represented the parametrically modulated 
activity by the risk level associated with each pump (i.e., the probability of explosion) and modeled all of the 
pump events at the onset time of pumping; and (2) Pumps,,,, the averaged activity across all pumps without 
parametric modulation by the individual average risk level (i.e., the mean pumps correspond to the probability of 
explosion) using the same onset time”. For a cash-out event, two regressors (similar to pumps) were defined: (1) 
Cash-Out arametrió and (2) Cash-out,,,. For an explosion event, two regressors similar to pumps were defined: (1) 
Explosionyarametrid and (2) Explosiongys. 

‘These six regressors were convolved with the gamma function to approximate the hemodynamic response 
of the brain. Head motion parameters in 6 directions were also included as covariates. The resultant z-statistic 
map (transferred from the t-value map from the GLM result) of PumpSparametric (i.€., the risk level) was trans- 
formed to the Talairach space. To test the activation of risk level, a group-level one-sample t-test was used and the 
responses were identified if they survived whole brain correction for family-wise errors at a cluster-level threshold 
ofp < 0.005 (cluster size 30 voxels, 810 mm?) and a voxel-level threshold of p < 0.0001. 


Risk level-related activation and relationship to task performance and severity of nicotine 
dependence. For the z-statistic map of risk level, the averaged z-value of every voxel within each signifi- 
cantly activated region related to risk level was calculated for each participant. Then, Pearson correlation analysis 
was used to assess the relationship between risk level-related brain activation and task performance as well as the 
severity of nicotine dependence. 

The same analysis was also applied to cash-out and explosion events to clarify to influence of outcome on risky 
decision making and the results were provided in the supplementary materials. 


RSFC and its relationship to task performance and severity of nicotine dependence. For the 
resting-state {MRI data, we regressed the motion data out of the time series and performed bandpass temporal 
filtering (0.01-0.0 Hz) on the residual signals. Then, to further reduce nuisance signals, the average white matter 
and CSF signals were regressed out. The white matter mask was determined from the high-resolution struc- 
tural image using the FAST segmentation program (Functional MRI of the Brain software library; www.fmrib. 
ox.ac.uk). Each participant’s CSF mask was manually drawn to fit the anatomical boundaries of a standardized 
Talairach atlas brain, which was transformed onto the image space of each participant and modified according to 
the cortical structures of the participant's brain by referencing the anatomical boundaries in the high-resolution 
structural image. These nuisance signals were used to account for fluctuations that were unlikely to be relevant to 
neuronal activity”. Then, data scrubbing was performed, and any volume with a framewise-dependent value for 
the temporal derivative of a time course’s root mean squared head motion variance exceeding 0.5 was excluded. 
The resultant resting-state fMRI data were then used for functional connectivity analysis. 

We defined the AI, dACC, and thalamus activated by risk level during BART as regions of interest (ROIs) in 
the RSFC analysis. The preprocessed resting-state fMRI time series were averaged within each ROI. Correlations 
between the averaged ROI time series were calculated and then transformed to Fisher z-values. Finally, the cor- 
relations between RSFCs and the mean pumps and the severity of nicotine dependence were calculated. We 
performed a FDR-correction for multiple comparisons (FDR < 0.05). 


Mediation analysis. To examine whether the RSFC of the dACC-thalamus pathway might mediate the effect 
of the FIND score on mean pumps, we performed a causal mediation analysis”*”>. First, the hypothesized medi- 
ator (the RSFC of the dACC-thalamus) was regressed on the independent variable (FTND score); second, the 
dependent variable (mean pumps) was regressed on the independent variable; and third, the dependent variable 
was regressed on both the independent variable and the hypothesized mediator in one equation. Coefficients for 
the equation were estimated, and mediation was considered to occur if the relationship between the independent 
and dependent variables could be partially or totally accounted for by the hypothesized mediator. 
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